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We extract the two-photon exchange (TPE) contributions to electron-proton elastic scattering 
using two parametrizations and compare the results to different phenomenological extractions and 
direct calculations of the TPE effects. We find that many of the extractions give similar results, 
and highlight the common assumptions and the impact of not including these assumptions. We 
provide a simple parametrization of the TPE contribution to the unpolarized cross section, along 
with an estimate of the fit uncertainties and the uncertainties associated with the model dependence 
of the extractions. We look at the contributions as extracted from various e-p elastic scattering 
observables, and make predictions for ratio of positron-proton to electron-proton elastic scattering 
cross sections. 
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I. INTRODUCTION 

Electron scattering is a powerful technique used to re- 
veal the underlying structure of the nucleon. This is 
because the electron is a point-like, structureless par- 
ticle. This makes it an ideal probe of the target, as 
the cross section represents only known couplings and 
the (unknown) structure of the target. In elastic elec- 
tron scattering, the incident electron scatters from a 
nucleon through the exchange of a virtual photon 7*, 
and the structure of the target nucleon appears through 
two Sachs [l| electromagnetic form factors, Ge p and 
Gmp, which measure the deviation of the scattering from 
that for a point-like spin- 1/2 target. As Q 2 , the four- 
momentum squared of the virtual photon, increases the 
scattering becomes more sensitive to the small scale in- 
ternal structure of the target. 

Utilizing electron scattering, there are two methods 
used to extract the proton's form factors. The first is the 
Rosenbluth or Longitudinal- Transverse (LT) separation 
method Q , which uses measurements of the unpolarized 
cross section, and the second is the polarization trans- 
fer or polarized target (PT) method [3], which requires 
measurement of the spin-dependent cross section. In the 
Rosenbluth separation method, the reduced cross section 
an is defined in the one-photon exchange (OPE) approx- 
imation as 



a R = G\ Ip {Q 2 ) + -G 2 Ep {Q 2 ), 



(1) 



where r = Q 2 /AM 2 , M p is the mass of the proton, and e 
is the virtual photon longitudinal polarization parameter, 
defined as e -1 = [l + 2(1 + r) tan 2 (^)] , where 9 e is the 
scattering angle of the electron. For a fixed Q 2 value, the 
reduced cross section an is measured at several e points, 
and a linear fit of an to e gives G 2 Mp {Q 2 ) as the intercept 
and G 2 Ep (Q 2 )/t as the slope. 



In the recoil polarization method, a beam of longitu- 
dinally polarized electrons scatters elastically from un- 
polarized proton target. The electrons transfer their po- 
larization to the unpolarized protons. By simultaneously 
measuring the transverse, Pt, and longitudinal, Pi, po- 
larization components of the recoil proton, one can de- 
termine the ratio h p Ge p /Gm p in the OPE 



G Ep _ Pt [E + E') 



Gm p 



Pi 2M P 



tan(y), 



(2) 



where E and E' are the initial and final energy of the in- 
cident electron, respectively. The ratio can be extracted 
in a similar fashion using polarized beams and targets 
by measuring the asymmetry for two different spin direc- 
tions @,0|. 

The two methods yield significantly and strikingly dif- 
ferent results in the region Q 2 > 1.0 (GeV/c) 2 , with val- 
ues of fj, p G e p I 'G m p differing almost by a factor of three at 
high Q 2 . In the LT separation method, the ratio shows 
approximate form factor scaling, h p Ge p /Gm p ~ 1, al- 
beit with large uncertainties at high Q 2 values. On the 
other hand, the recoil polarization method yields a ra- 
tio that decreases with increasing Q 2 . The ratio is well 
parametrized Q as 



H p Ge p 
Gm p 



H P R = 1 - 0.13(Q 2 - 0.04). 



(3) 



To reconcile the ratios, several studies suggested that 
missing higher order radiative corrections to the electron- 
proton elastic scattering cross sections can explain the 
discrepancy. In particular, studies focused on the role of 
the two-photon exchange (TPE) effect in resolving the 
discrepancy. The TPE effect was studied extensively 
both theoretically [9|4l4| and phenomenologically [la - 
[23j . An extensive review of the role of the TPE effect 
in electron-proton scattering can be found in 0, [H[ . 
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II. 



TWO-PHOTON-EXCHANGE AND THE e ± p 
CROSS SECTION SCATTERING RATIO 



The interference of the OPE and TPE amplitudes rep- 
resents the leading TPE correction to the electron-proton 
elastic scattering cross section. However, the recoil po- 
larization data were confirmed experimentally to be es- 
sentially independent of e [26| . We account for the TPE 
contribution to gr simply by adding the real function 
F(Q 2 ,e) to the Born reduced cross section: 



<TR = G 2 Mp (l+ e -R 2 ) + F(Q 2 ,e) 



(4) 



where R = Ge p /Gm p is the recoil polarization ratio Eq. 
(|3]). The TPE correction changes sign for electron-proton 
and positron-proton scattering, so the ratio R e + e ~ (Q 2 ,e) 
defined as 



R, 



\Q\e) 



<j(e + p 



a(e p 



-pV 



(5) 



allows for direct measurements of the TPE contribution 
to the cross section. The ratio R e + e ~ can be expressed 
as R e + e - »1 + 45ft(A 27 )/Ai 7 , with and A 2l being 
the OPE and TPE amplitudes [27j, respectively. Here 5R 
stands for the real part. On the other hand, the modifi- 
cation to the electron cross section is 1 — 25R(A 27 )/A l7 . 
Clearly, any change in the electron cross section will have 
almost twice the change in the ratio R e + e - but with the 
opposite sign. 

Recently, several theoretical studies estimated the 
function F(Q 2 , e) fijl-Hij. Experimentally, studies focused 
on searching for nonlinearities in <jr [16l4l8l l28j and phe- 
nomenologically by reanalyzing the experimental data 
usin g a proposed p arametrization of the TPE contribu- 
tion [LQ, \WH23^ . We present a summary of the results 
which we will examine in detail in this work. 

Guichon and Vanderhaeghen [l5j expressed the 
hadronic vertex function in terms of three independent 
complex amplitudes (generalized form factors) which de- 
pend on both Q 2 and e: Ge p (£,Q 2 ), Gm p (^,Q 2 ), and 
i*3(e, Q 2 ). These generalized form factors can be broken 
into the usual Born (OPE) and the TPE contributions 

as: Ge p ,m p (£,Q 2 ) — Ge p ,m p {Q 2 ) + AGe p ,m p {s,Q 2 ) 
with Y 2l {v, Q 2 ) defined as ^( M ajg 3 fl/ | ) , where v — 

M£y/(1 + e)/(l - e)y/T(l + t). The reduced cross sec- 
tion is expressed in terms of these amplitudes as 



<TR= \G 



M P \ 



1 £ \G Ep \ 2 



T \Gm p 



-2e 1- 



\G 



Ep\ 



t\Gm p \ 



Y 
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(6) 



They demonstrated that small TPE contributions could 
significantly modify the extraction of the form factor 
ratio from the Rosenbluth separation while having rel- 
atively little impact on the polarization transfer data. 
They also provided an estimate of the TPE amplitudes, 
based on the assumption that the entire contribution 



comes from the Y 2l (v, Q 2 ) term, parametrized so as to 
yield a correction to or that is proportional to e. This 
yields a small reduction in Ge p as extracted from the 
polarization measurement, and a significant reduction in 
the Rosenbluth extraction of Ge p , such that it becomes 
consistent with the corrected polarization results. In this 
work, Gm p is unaffected by the TPE contributions. 

Based on the framework of [HI], Arrington [lj| per- 
formed a global analysis where he extracted the TPE 
amplitudes AGe p ,m p and Y"2 7 • He assumed that the am- 
plitudes were e-independent and took AGe p — 0. Values 
for l2 7 (Q 2 ) were extracted from the difference between 
polarization and Rosenbluth measurements, taking into 
account the uncertainties in both data sets. Based on 
the high-e constraints from the comparison of positron 
and electron scattering [29| the amplitude AGm p was 
determined by requiring that its contribution to <jr at 
e = 1 cancelled the contribution of Y2 7 . The extracted 
TPE amplitudes and their estimated uncertainties are 
then parametrized as a function of Q 2 , and used to ap- 
ply TPE corrections to the form factors obtained from a 
global Rosenbluth analysis [3(| and the new recoil polar- 
ization data. Throughout this text, we will refer to the fit 
to the uncorrected Rosenbluth form factors obtained in 
Ref. [13] as the Arrington cfr Fit, and to those corrected 
using the extracted TPE amplitudes as the Arrington Y 2l 
Fit 111. 

In both of these analysis, there is not enough informa- 
tion to directly determine the amplitudes, so assumptions 
have to be made about the relative importance and the 
e-dependence of the three TPE amplitudes. Common 
to these analyses are the assumption that the correction 
is close to linear in e, as no non-linearities have been 
observed [H-[ll HI]. If one also neglects the TPE cor- 
rection to the polarization data, which is significantly 
smaller at high Q 2 , then it is not necessary to work 
in terms of the polarization amplitude; one can simply 
parametrize the TPE contributions to the reduced cross 
section, taking a linear (or nearly linear) e dependence. 

Alberico et al. [3l| performed a global extraction of 
the proton form factors and TPE contributions, based 
on one of the parametrizations of the TPE contributions 
from Chen et al. uM: 



OR = G 2 



Alp 



1 + -R 2 

T 



+ A{Q 2 )y + B{Q 2 )y\ (7) 



where y = ^(1 - e)/(l + e), A(Q 2 )=aG 2 D (Q 2 ), 
B(Q 2 )—PG 2 D {Q 2 ) and Gd(Q 2 ) is the dipole parameteri- 
zation: Gd(Q ) = [l + Q 2 /(0.71(GeV/c) 2 )]- 2 . Note that 
Chen et a/.performed such an extraction, but analyzed 
only the data of Ref. [321 ]. Alberico et al. extract the form 
factors and the parameters a and j3 in two different ways. 
In the first fit, referred to throughout this text as ABGG 
Fit 1, the polarization transfer data was used to fix the 
ratio R(Q 2 ) = fipG Ep /G Mp = 1.022 - 0.13Q 2 and G Mp 
and the TPE contribution from Eq. were determined 
in a global fit to world's cross section data. Gm p was 
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parametrized using the functional form from Kelly [3J 
G Mp (Q 2 ) 1 + ELi < k r k 



(8) 



with n = 1, yielding four parameters for Gm p and the two 
TPE parameters a and /3. In their second fit (ABGG 
Fit 2), the cross section and polarization transfer data 
were fit simultaneously, using 4 parameter fits to Ge p 
and Gmp along with the two TPE parameters. 

Qattan and Alsaad [21| proposed a different empirical 
parametrization for the function F(Q 2 ,e). The function 
F(Q 2 , e) was double Taylor series expanded as a polyno- 
mial of order n keeping only terms linear in e but without 
constraining the TPE amplitudes by enforcing the Regge 
limit. In their final result, labeled fit III, the reduced 
cross section ctr was parametrized as 



a R — G Mp ( 1 + -R 



ef(Q 2 ), 



(9) 



where for the extraction from a fixed-Q 2 data set, f(Q 2 ) 
is just a constant. Throughout this text, Eq. (j9|) will be 
referred to as the QA parametrization. 

As noted in some previous extractions [ltl [2(| [2^, [2t| , 
small angle comparisons of positron to electron scatter- 
ing set tight limits on the TPE effects and TPE calcu- 
lations [J [HI, HH Hfl show that the contribution should 
be zero at e = 1. Therefore, we perform an updated 
version of the fit from Ref. [21| using a parametriza- 
tion which maintains the linear correction in e but yields 
F(Q 2 ,e) = at e = 1. Based on the parametrization 
from Borisyuk and Kobushkin [22| . we take the following 
form 



<?R 



r 2 

^Mp 



1 



-R' 

T 



2o(l - e)G 2 Mp , 



(10) 



for our updated extraction. Throughout this text, 
Eq. (fTU| will be referred to as the BK parametrization. 
Because the recoil polarization ratio Ge p /Gmp was ex- 
perimentally confirmed to be independent of e (26|, the 
ratio R = Ge p /Gmp was fixed to be that of the recoil 
polarization ratio or Eq. ©, as in the analysis of Ref. [2l[ 



III. RESULTS AND DISCUSSION 



In this article we do the following: 
(1) We extend the analysis of Ref. 



211 ] . extracting the 



TPE contributions to or using both the QA and BK 
parametrizations for the TPE amplitude. We also include 
additional data sets [33 - l36j , that were not included in the 
original analysis. For the data from Qattan et al. [37j 
we use the published results, including the e-dependent 
systematic uncertainties, yielding somewhat modified re- 
sults than in the previous analysis (2lj . 

(2) We take the proton form factors and TPE ampli- 
tudes extracted using the QA and BK parametrizations 
and compare them to those obtained using the extracted 



TPE corrections Refs. [19|, |3fl |3l| and the TPE calcula- 
tion of Ref. [38| ■ We also use our results and the previous 
parametrizations to determine the ratio R e + e -. We use 
these results to examine the Q 2 dependence of the ex- 
tracted TPE corrections, the consistency of the different 
extractions, and discuss the impact of the different as- 
sumptions used in these analyses. 

(3) We compare the TPE contributions as extracted 
using these different approaches to more recent extrac- 
tions [2(], [22| which attempt to extract the TPE ampli- 
tudes using additional constraints that come from a re- 
cent measurement of the e dependence of Ge p /Gm p at 
Q 2 = 2.50 (GcV/c) 2 H|. ' 



A. Form Factors and the TPE Amplitudes 

We begin with the extraction of the form factors and 
TPE contributions based on the QA parametrization, 
Eq. dp. The results of this fit for the data sets of 
Refe. gUMMIMSl are given in Table H The extrac- 
tion follows the procedure of fit III from Ref. [2l(, with 
the fit to the polarization transfer data, Eq. used to 
constrain R, leaving Gm p and / as the two fit parame- 
ters. We also extract the TPE contributions using the 
BK parametrization Eq. (ITU|) . Table HT1 lists the results 
of the BK fit. 

Note that for these older data sets, we take the cross 
sections used in the analysis of Refs. [3(| 53, where miss- 
ing higher older radiative corrections terms such as the 
Schwinger term and additional vacuum polarization con- 
tributions from muon and quark loops have been applied. 
See Ref. 42j for tabulated electron-proton elastic scatter- 
ing cross sections data. 

Figure [1] shows the values of (Ge p /Gd) 2 and 
(G m p I ' UpG d) 2 obtained using the QA and BK 
parametrizations (shown as hollow and solid circles, re- 
spectively). The solid curves are fits to the extracted 
form factors from previous phenomenological analyses: 
Ref. [H| (ABGG Fit I and ABGG Fit II), and Ref. [H 
(Arrington Y27 Fit), and the form factors extracted based 
on a calculation of the TPE contributions from Ref. 38 1 
(AMT Fit). The dotted line is the global Rosenbluth ex- 
traction from Ref. [3(| (Arrington or Fit), which does 
not include TPE corrections. 

The QA results for Gup are typically 2-10% below 
the BK values. This simply reflects the fact that the QA 
parametrization yields no TPE contribution at e = 0, 
and thus the (linear) extrapolation to the anchor 'point 
at e = is unchanged by the application of the extracted 
TPE contributions, while the extrapolation is modified 
by the TPE contributions for the BK parametrization. 
Note that because they use the same PT parametriza- 
tion for UpGEplGMp, the Ge p values are also lower by 
the same amount for the QA extraction. In addition, 
the uncertainties are much larger in the QA extraction 
for some kinematics. Most experiments tend to have a 



large part of their data at large e values, so there is a 
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TABLE I: The form factors and TPE parameters obtained 
using the QA parametrization, Eq. @, as a function of Q 2 
(given in units of (GeV/c) 2 ). The TPE correction is ex- 
tracted based on the assumption that it fully resolves the 
difference between /j. p Ge p /Gm p as extracted from the given 
Rosenbluth extraction and the value from the polarization 
transfer parametrization of Eq. ©. The overall normaliza- 
tion uncertainty for each data set (typically 1.5-3%) is not 
taken into account in this extraction. See text for complete 
details. 

Q 2 {Gm p /li p G d ) 2 {Ge p /G d ) 2 f(Q 2 ) x 100 x l 
Andivahis [32] (taken from Rcf. [21]) 
1.75 1.106±0.014 0.669±0.009 0.32±0.14 0.30 
0.514±0.006 0.07±0.04 0.53 
0.371±0.006 0.05±0.02 0.16 
0.251±0.004 0.02±0.01 0.51 



2.50 
3.25 
4.00 
5.00 



1.111±0.013 
1.092±0.018 
1.068±0.017 
1.029±0.017 



0.130±0.002 0.01±0.005 0.93 



Walker [39] (taken from Ref. [21J) 

1.00 1.024±0.055 0.785±0.042 2.40±1.97 0.72 

2.00 1.020±0.037 0.563±0.020 0.67±0.18 0.65 

2.50 1.040±0.034 0.483±0.016 0.28±0.08 0.98 

3.00 1.010±0.045 0.378±0.017 0.18±0.06 0.19 



Christy [40J (taken from Ref. [21]) 



0.65 
0.90 
2.20 
2.75 
3.75 
4.25 
5.25 



0.916±0.061 
1.047±0.036 
1.110±0.030 
1.115±0.021 
1.087±0.030 
1.023±0.024 
1.015±0.065 



0.777±0.052 
0.826±0.028 
0.574±0.016 
0.468±0.009 
0.291±0.008 
0.210±0.005 
0.106±0.007 



13.07±6.33 0.02 
1.58±1.81 1.37 
0.13±0.12 
0.07±0.05 
0.03±0.02 
0.05±0.01 
0.02±0.01 



1.07 
0.04 
1.15 
0.55 
0.78 



Qattan [37] 
2.64 1.108±0.007 
3.20 1.098±0.008 
4.10 1.064±0.010 



0.486±0.003 
0.381±0.003 
0.237±0.002 



6.39±1.10 0.35 
8.13±1.32 0.54 
10.80±1.78 0.14 



Bartel [34] 
1.169 1.045±0.029 
1.75 1.103±0.029 



0.761±0.021 
0.667±0.018 



1.63±0.89 0.26 
-0.03±0.30 0.01 



Litt [35] 

1.50 0.970±0.247 
2.00 0.961±0.230 
2.50 1.031±0.062 
3.75 0.960±0.081 



0.637±0.162 
0.534±0.128 
0.477±0.029 
0.257±0.022 



1.86±2.67 0.01 

0.86±1.09 1.10 

0.23±0.15 0.65 

0.10±0.05 0.32 



Berger [36J 

0.389 0.970±0.037 

0.584 0.970±0.017 

0.779 1.008±0.025 

0.973 1.005±0.032 

1.168 1.047±0.048 

1.363 1.076±0.046 

1.558 1.063±0.057 

1.752 1.129±0.077 



0.884±0.034 
0.837±0.015 
0.823±0.020 
0.776±0.025 
0.763±0.035 
0.738±0.031 
0.685±0.036 
0.682±0.047 



-0.9±12.8 
5.27±2.78 
1.82±2.22 
3.37±1.65 
2.35±1.94 
0.29±1.22 
1.71±1.48 
0.07±1.21 



0.48 
0.38 
0.55 
0.81 
0.29 
0.58 
1.12 
0.75 



TABLE II: The form factors and TPE parameters obtained 
using the BK parametrization, Eq. (jTUJ) , as a function of Q 2 
(given in units of (GeV/c) 2 ). The TPE correction is ex- 
tracted based on the assumption that it fully resolves the 
difference between ll p Ge p /Gm p as extracted from the given 
Rosenbluth extraction and the value from the polarization 
transfer parametrization of Eq. The overall normaliza- 

tion uncertainty for each data set (typically 1.5-3%) is not 
taken into account in this extraction. See text for complete 
details. 
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FIG. 1: (Color online) (Ge p /Gd)^ [topland (G Mp / H P G D ? 
[bottom] as obtained from Refs. [3^. l34 - l37l . [39. 40] using the 
QA and BK parametrizations. Also shown are the ABGG I 
and ABGG II fits 0, the AMT fit [aj|, the Arrington Y 2l 
Fit , and the Arrington an Fit [3(| which does not include 
any TPE contributions or polarization data. 

significant extrapolation to e — in the QA extraction, 
yielding a larger uncertainty in the overall normalization. 
For the BK parametrization, the normalization is fixed 
at the s = 1 value. This reduction in the uncertain- 
ties demonstrates one of the important strengths of the 
polarization transfer measurements, as they significantly 
reduce the uncertainty associated with the extraction to 
e = 0. While the uncertainty on the PT data is small, 
neglecting this uncertainty in the fit (as we and some 
other extractions do), yields a small underestimate of 
the uncertainty. The uncertainties associated with this, 
and other assumptions that go into the phenomenological 
extractions, will be discussed in Sec. IIII B[ 

In the extraction of (Ge p /Gd) 2 , there is a large dif- 
ference between the TPE-uncorrected result (Arrington 
<tr) and the extractions that apply corrections for TPE 
contributions. This is because the LT extraction of Ge p 
at large Q 2 is extremely sensitive to angular-dependent 
TPE corrections, since Ge p enters into the cross section 
as a small, e-dependent term. Nearly all of the other 



fits extractions are in excellent agreement, but this is 
simply because they all explain the difference between 
the LT and PT extractions of Ge p /Gm p in terms of 
a nearly-linear correction to the reduced cross section. 
Therefore, these analyses will, by construction, ensure 
that the final result for Ge p /Gm p will be consistent with 
the polarization data, thus yielding nearly identical val- 
ues for Ge p except for small differences in the extracted 
values of {Gm p I H p Gd) 2 , most clearly seen in the QA 
extraction. The quality of the agreement shows that 
the small non-linearities in the ABGG analysis and the 
choice of cross section data sets and parametrization for 
the polarization transfer measurements has little impact. 
The AMT analysis applies TPE contributions based on 
a hadronic calculation (Toj |. with a small additional con- 
tribution designed to more fully resolve the discrepancy 
for Q 2 > 2 (GeV/c) 2 , where the calculated TPE do not 
bring the LT and PT results into perfect agreement. Be- 
cause of this additional contribution, and the neglect of 
TPE contributions to the polarization data, this analysis 
must also yield results consistent with the other analy- 
ses at large Q 2 . The only analysis that yields different 
results is the Arrington Y27 fit, where the extraction of 
all three amplitudes allows for a contribution of the TPE 
contribution to the polarization transfer data, yielding a 
noticeable downward shift in the value of h p Ge p /Gm p 
and thus (Ge p /Gd) 2 - However, the size and even the 
sign of this correction depend on the assumptions made 
in trying to separate the three TPE amplitudes. Thus, 
the deviation from the other fits is at best an indication 
of the possible uncertainty in these extractions. This will 
be addressed further when we compare to other extrac- 
tions that include TPE contributions to the polarization 
data in Sec. ITITbI 

For the extraction of (Gm p / H p Gd) 2 , the BK results 
are in good agreement with previous extractions that 
include TPE, while thc QA results are more consistent 
with the Arrington ctr fit which neglects TPE. The BK 
and ABGG other extractions assume no TPE contribu- 
tion to the polarization data and an approximately lin- 
ear correction to the reduced cross section, as discussed 
above. Because they require very similar corrections to 
the e dependence of <jr to resolve the discrepancy in 
H p Ge p /Gm p measurements, they have very similar cor- 
rections to (Gm p I H p Gd) 2 , which simply relates to the 
low-e value of the correction. The Arrington Y27 fit is 
the only extraction in which the polarization transfer re- 
sults for Ge p /Gm p are modified by TPE contributions. 
While the overall size of the TPE correction required to 
resolve thc LT-PT discrepancy is somewhat larger in this 
case, the absolute correction to the cross section is small, 
yielding very similar values for (Gm p / ^ p Gd) 2 ■ 

Figure [2] shows the fit parameter a(Q 2 ) as a function of 
Q 2 for all data sets. The parameter a(Q 2 ) is at the few 
percent level, and for the most part, increases in magni- 
tude with increasing Q 2 . Note that the TPE contribution 
to the cross section is F(Q 2 ,e) = 2a(l~e)G 2 Mp , and thus 
at high Q 2 , the fractional slope introduced by the TPE 
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correction is 2a. The extracted TPE contributions from 
the different data sets are relatively consistent, and show 
a slow increase in the TPE contribution as Q 2 increases. 



0.1 



0.05 



2 

~9 



0.05 



-0.1 



Bartel 
□ Litt 
Berger 

- BMT Correction 

- AMT Correction 




U - A n) 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 



Q 2 [(GeV/c) 2 ] 



FIG. 



2: (Color online) The fit parameter a(Q 2 ) as obtained 
using the BK parametrization from the data of Refs. [3^ . l34r - 
H3i [H, 51}. The dashed curves correspond to fits to the Q 2 
dependence using different parametrizations (see text) and 
the solid lines correspond to the values of a(Q 2 ) determined 
by fitting to the TPE calculations of Ref. [I(J, and the TPE 
correction of Ref. i38j which adds and additional phenomeno- 
logical TPE contribution at higher Q 2 . 

In an attempt to parametrize the Q 2 dependence of the 
parameter a(Q 2 ), several different functional forms were 
tried and are plotted as dashed lines in Fig. [2] These fits 
all give reasonable reduced x 2 values, 1.04 < xt < 1-10. 
The lowest xt value was obtained using the form a(Q 2 ) — 
a \fQ 2 with a=— 0.0191±0.0014 (magenta dashed line in 
Fig. E}. The full spread of the fits is below 0.005 for 
1<Q 2 <4 (GeV/c) 2 , although at lower Q 2 values, none 
of the extractions are precise and the behavior is largely 
unconstrained below Q 2 w 1.5 (GeV/c) 2 . For our global 
fit, we take 

a(Q 2 ) = -0.0191 y^Q 2 ± 0.0014a/Q2 ± 0.003 (11) 

where the first uncertainty is the fit uncertainty in a, 
and the second is the systematic error band included to 
account for the model dependence of the fit for 1.5 < 
Q 2 < 4 (GeV/c) 2 . 

Figure [2] also shows curves for a[Q 2 ) based on the 
hadronic TPE correction of Ref. along with the 

version used in Ref. [38j which includes a small addi- 
tional contribution at high Q 2 values. The version used 
in the global analysis of the form factors 38] is in good 
agreement with our fit for Q 2 J> 2 (GeV/c) 2 , where the 
data provides significant constraints on the TPE contri- 
butions. Note that the calculated TPE corrections, as 
well as those of Refs. [U HJ, show a change of sign in 
the TPE effects for Q 2 < 0.5 (GeV/c) 2 . This is not seen 
in the data, although the constraints at low Q 2 are in- 
sufficient to make strong conclusions in this region. 



The fact that the low Q 2 behavior is inconsistent with 
the calculations is not entirely surprising. Our extrac- 
tion, like most similar phenomenological analyses, as- 
sumes that TPE contributions are significant for the 
cross section measurements but negligible for polariza- 
tion data. However, both the hadronic [10J and par- 
tonic [45] calculations suggest that the TPE contribu- 
tions are at the few percent level for both observables. 
The main difference is that the impact of the correction 
on the extracted form factors is amplified for G e p at high 
Q 2 , where the form factor is extracted from a small angu- 
lar dependence in the cross section which can be notice- 
ably modified by a few percent TPE contribution. At low 
Q 2 , the TPE contribution is not amplified in the extrac- 
tion of the form factors from the cross section data, so 
neglecting the contributions to the polarization transfer 
data will not be reliable. In addition, because the dif- 
ference between h p Ge p /Gm p from Rosenbluth and po- 
larization transfer measurements becomes small at low 
Q 2 , neglecting the uncertainty in the polarization trans- 
fer extraction of h p Ge p /Gm p leads to a significant un- 
derestimate of the uncertainties below 1 (GeV/c) 2 . 




FIG. 3: (Color online) The Ratio F(Q 2 , e)/G 2 Mp as a function 
of e at Q 2 = 1.75, 2.50, 3.25, 4.0, and 5.0 (GeV/c) 2 from 
Ref. [13] determined using the QA parametrization (magenta 
line), the BK parametrization (cyan line), the ABGG Fit I 
(black line) and the ABGG Fit II (red line), and the Arrington 
Y 2j Fit (blue line). 

In Fig. [31 we compare the e dependence of the TPE 
corrections based on the different extractions based on 
the data from Ref. [3_2j. In addition to the extractions 
from the QA and BK parametrizations, we show the 
ABGG Fit I (black line), ABGG Fit II (red line), the 
Arrington Y27 Fit (blue line). Most of the extractions 
yield similar slopes, with the QA extraction differing 
from the others in that the TPE contribution goes to 
zero at s = 0. The Arrington Y27 parametrization has a 



7 



somewhat larger slope, because of the inclusion of TPE 
contributions which reduce R as measured in polariza- 
tion experiments, thus necessitating a greater decrease 
to R as extracted in a Rosenbluth separation. The TPE 
corrections increase slowly with Q 2 except for the ABGG 
fits which show essentially no Q 2 dependence. In these 
extractions, the Q 2 dependence of the TPE contribution 
is taken to go as G 2 J (Q 2 ), and thus they are a nearly con- 
stant fractional correction to the cross section at large Q 2 
values, where Gm ~ l^pGn dominates the cross section. 



B. The R e + e - Ratio 

The function F(Q 2 ,e) which represents the interfer- 
ence of the OPE and TPE amplitudes, changes sign de- 
pending on the charge of the projectile, yielding an ampli- 
fied signal when taking the ratio of electron and positron 
scattering. The ratio R e + e - {Q 2 7 e), defined by Eq. (|5|), 
is determined simply by changing the sign in front of the 
TPE amplitudes. We determined the ratio R e + e - using 
the form factors and the TPE amplitude extracted us- 
ing the QA and BK parametrizations. These ratios will 
then be compared to those obtained from other analy- 
ses [H US Hi, ill, and calculations HI El ■ We also 
use our extracted TPE contribution to make predictions 
for new and recently completed measurements [46|-|50(. 
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FIG. 4: (Color online) R e + e - as a function of e at Q 2 = 1.75, 
2.50, 3.25, 4.0, and 5.0 (GeV/c) 2 from Ref. [H] determined 
using the QA parametrization (solid black circles), and the 
BK parametrization (open red squares). The solid black and 
red lines through the data are calculated using Eqs. (|12[) and 
(|13J1 . respectively, while the dashed red line is the global fit 
from the BK extraction Eq. Also shown are results of 

the Arrington Y 2l Fit (solid blue line), the ABGG Fit I (black 
dashed line) and the ABGG Fit II (black dotted line). 



e extracted from the data of Ref. [32[ using the 
QA parametrization (solid black circle) and the BK 
parametrization (open red squares), along with the pre- 
vious extractions. The ratio R e + e - as determined using 
the QA parametrization is always equal to or less than 
unity, while the other parametrizations force R e + e - = 
at e = 1, yielding a ratio that is always equal to or larger 
than one. Because measurements of R e + e ~ put signifi- 
cant constraints on deviations from unity at e = 1, we 
exclude the QA parametrization from further compar- 
isons. As seen in Figs. [3] and 21 the e dependence is very 
similar to that of the BK fit, except for the overall offset. 

For a better comparison between the two parametriza- 
tions, we fit the ratio extracted using the QA 
parametrization to the form 



Rp 



\Q 2 ,e) = l + B{Q 2 )e 1 



(12) 



and that extracted using the BK parametrization to 



R e+e -(Q 2 ,e) = l-B{Q 2 )(\-e), (13) 



with B{Q 2 ) being the parameter of the fit and represents 
the slope. The values of B(Q 2 ) from the BK fit are given 
in Table UTT1 and shown as a function of Q 2 for each exper- 
iment in Fig. [5j The slope B(Q 2 ) is negative and grows 
in magnitude with increasing Q 2 value almost for most 
data sets. We also compare the results to a reference 
curve corresponding to our global fit Eq. (1TTT) . 
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FIG. 



5: (Color online) The slope B(Q ) as a function of Q 
as determined using the BK parametrization. The solid lines 
connecting the data points are to guide the eye, while the 
black dashed curve is our global best fit. 

In Figs. S] and El we compare previous TPE extractions 
to our results for the kinematics of the two highest preci- 
sion data sets [H, The ratio R e + e - as predicted by 
ABGG Fit 1 and ABGG Fit 2 are nearly identical, but 
Figure @] shows the ratio R e + e - as a function of smaller than the Arrington Y27 Fit. For the Andivahis 
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TABLE III: The values of the slope B(Q 2 ) as obtained using 
the QA and BK parametrizations by fitting to Eqs. (|12[1 and 
respectively. 
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data, our extraction is in generally good agreement with 
the ABGG fits, but the Qattan data yield larger TPE 
corrections, especially at their higher Q 2 values. 

At present, precise measurements of R e + e - are limited 
to relatively low Q 2 or large e, making it difficult to di- 
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FIG. 6: (Color online) R e + e - as a function of e as extracted 
from Ref. [37| using the BK parametrization (square points 
and solid red line). The red dashed lines show the global fit 
of Eq. pip . Also shown are the results of the Arrington Y27 
Fit (solid blue line), the ABGG Fit I (black dashed line) and 
the ABGG Fit II (black dotted line). 



rectly compare the data to estimates of the e dependence 
of TPE at high Q 2 extracted from comparisons of the 
Roscnbluth and polarization measurements. As stated in 
Ref. [29| . if TPE corrections are responsible for the dis- 
crepancy between the Rosenbluth and recoil-polarization 
data, a 5-8%, linear or quasi-linear, e-dependent correc- 
tion to the electron cross section is required to resolve the 
discrepancy at high Q 2 . This implies that R e + e - at high 
Q 2 should have a 10-16% quasi-linear e dependence, with 
R e + e - decreasing as e increases. The existing R e + e - data 
provide some evidence for such a e dependence, but the 
significance is only 3ct and limited to low Q 2 data, where 
the e-dependence is observed to be (5.7±1.8)% [29]. This 
corresponds to B(Q 2 ) = 0.057, but the extraction uses 
only data below Q 2 = 2 (GeV/c) 2 , corresponding to an 
average Q 2 values of approximately 0.5 (GeV/c) 2 , mak- 
ing the result consistent with the low Q 2 extractions pre- 
sented here. Relatively precise data do exist at large e 
and moderate Q 2 values, suggesting that R e + e - ~ 1 for 
e — > 1, as assumed in most of the extractions. 

Finally, we make an additional comparison of the e de- 
pendence at Q 2 = 2.64 (GeV/c) 2 , including two recent 
extractions which take advantage of new measurements 
of the e dependence of the recoil polarization data [2(| . 
These are similar to the earlier attempts to separate the 
individual TPE amplitudes [Hill, 

but because of the 
new data on the e dependence of the polarization mea- 
surement of h p Gep/Gmp and the polarization compo- 
nent Pi, fewer assumptions have to be made as to the 
contribution of the different amplitudes. 

Guttmann et al. [201 ] used the measurements of the £ 



9 





Our Fit Arlington Y 2y Fit 


1.15 


v ABGG I Fit 




- N. ABGG II Fit 




\_ Guttmann Fit 




Borisyuk and Kobushkin Fit 




NJ\ n. Blundenetal. 


1.1 


■= x. n. Arrington et al. 


+ 


**X >v Coulomb Distortion 






X 




1.05 






1,1,1,1 /"^i. 



h 0.2 0.4 0.6 0.8 1 



e 

FIG. 7: (Color online) Ratio R e + e - as a function of e at 
Q 2 = 2.64 (GeV/c) 2 . In addition to the curves shown 
in Fie. [6l we also show the extractions from Guttmann et 
al. [20] and Borisyuk and Kobushkin [22], and calculations 
of Coulomb distortion [3] and hadronic TPE contributions 
from Blunden et al. and with the addition of a small ad- 
ditional contribution at higher Q 2 values [3g]. Note that the 
extraction of Ref. [22| was performed at Q 2 = 2.50 (GeV/c) 2 . 
The uncertainties in our fit Eq. correspond to an uncer- 
tainty of ±0.015 at e=0. 

dependence of Pi and Pt/Pi J2al along with the e de- 
pendence of the cross section [371 ] to constrain the three 
TPE amplitudes Ym, Ye, and Y 3 . By writing the ob- 
servables in terms of Born values and the TPE ampli- 
tudes and keeping only the dominant TPE corrections 
(i.e. neglecting those that are suppressed by kinematic 
or other factors), they observe that the corrections to 
Pi/PP° rn are largely determined by Y3 and the correc- 
tions to Pt I Pi (which shows no e-dependence) come from 
a sum of I3 and Ye- Thus, the e dependence of P\ is 
used to determine Y3 , taking two different e dependences 
based on cither constraints from a pQCD calculation [51] 
or a modified parameterization. The lack of TPE con- 
tribution to Pt/Pi is then used to constrain constrains 
Ye- This allows Ym to be extracted based on the dif- 
ference between the predicted and measured values of 
the reduced cross section as a function of e. This ac- 
counts for all of the uncertainty in the extraction of the 
PT result for fj, p G e p / 'G m p based on the assumption that 
there is no £-dependent TPE correction, but does not ac- 
count for the fact that the data are still consistent with 
a e dependence of a few percent over the full e range. 
Their extracted amplitudes are used to determine R e + e - 
at Q 2 = 2.64 (GeV/c) 2 , which reaches a maximum value 
of 1.11 ± 0.016, is compared to the other extractions in 

Fig. m 

The analysis of Borisyuk and Kobushkin [22[ takes a 
similar approach, although they use a different linear 



combination of amplitudes than Ref. 1151 . Again, the e 
dependence of P t /Pi is taken to be zero [26| , and the cor- 
rection to the cross section is taken to be linear [TtJ ■ Data 
from available electron-proton scattering cross sections in 
the range of 2.20 < Q 2 < 2.80 (GeV/c) 2 were interpo- 
lated to Q 2 = 2.5 (GeV/c) 2 to extract the amplitudes 
at a fixed Q 2 value. This yields an extraction in terms 
of a single amplitude, and using the parametrization of 
Eq. CEO]), they obtain a value of a = -0.0250 ± 0.0035, 
corresponding to a peak R e + e - value of 1.11 ± 0.016. 

Note that the values in our result and both of these 
extractions of the TPE amplitudes are very similar, as is 
our uncertainty compared to the Guttmann et aLresult. 
While the details of the assumptions made in separating 
the TPE amplitudes affect the prediction they make for 
polarization observables, they both assume that there is 
no TPE contribution to the polarization transfer ratio, 
and thus are extracting amplitudes designed to induce 
a linear correction to the reduced cross section to make 
it consistent with the PT data. So for the observables 
related to the form factor extraction, this is the same 
assumption made in all of the phenomenological extrac- 
tions except for the Arrington Y27 fit. Similarly, these 
extractions typically neglect some of the uncertainties 
associated with the TPE contributions to the PT mea- 
surements. All of them ignore the uncertainty in the ex- 
tracted value of h p Ge p /Gmp or the impact of a possible 
e dependence, and most of the phenomenological extrac- 
tions ignore both. The Arrington Y27 fit includes the 
uncertainty in the extraction of n P G e p / 'G m p from PT 
data, but not the uncertainty associated with the TPE. 
One can take the difference between the Arrington Y27 fit 
and the other extractions of Ge p , Fig. [TJ as an estimate 
of the latter uncertainty, but the TPE impact on the PT 
data in this analysis is rather large, so this difference is 
probably better treated as an upper limit in the uncer- 
tainty, at least for Q 2 up to 3-4 (GeV/c) 2 . A shift in the 
PT values of fi p G Ep /G Mp of 0.02 at Q 2 w 1 (GeV/c) 2 
and 0.05 for Q 2 « 5 (GeV/c) 2 would yield a modified 
TPE correction, changing the extracted value of a(Q 2 ) 
by 0.003-0.005 (and the low e value of R e+e - in Fig.[7]by 
0.01-0.02), which is typically comparable to the total un- 
certainty we quote in the fit Eq. (|11|) and even larger for 
the lower Q 2 values of our fit. This shift in [i p Ge p IGm p 
is at the level of the uncertainties in the polarization ex- 
traction. It is also comparable to the level at which the £ 
dependence of the TPE contributions to the polarization 
data are constrained. Thus, the uncertainty neglected 
when the uncertainties and/or TPE contributions to the 
PT data are neglected are comparable to the total un- 
certainties obtained in these extractions. 

There has been a recent push to make new measure- 
ments of the ratio R e + e - , focusing on small e, where 
TPE contributions appear to be largest. The first is 
the VEPP-3 experiment [46T4481 ] . where the internal tar- 
get at the VEPP-3 electron-positron storage ring at 
Novosibirsk was used to extract the ratio R e + e - at 
Q 2 = 1.60 (GeV/c) 2 and e « 0.4. A raw R e+e - ra- 
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tio of 1.056±0.011 was obtained, which must be reduced 
by the charge-dependent bremsstrahlung correction, es- 
timated to be ~3%. The Blunden et al. calculation 
predicts R e + e - — 1.036 (1.043 with the additional con- 
tribution included in Ref. |38|). Our global fit predicts 
R e + e - — 1.060 ± 0.009, but the Q 2 value is low enough 
that the extraction is not expected to be very reliable, 
and the uncertainty associated with the model depen- 
dence of the extraction is larger than the quoted fit un- 
certainty. 

The second experiment is Jefferson Lab experiment 
E07-005 where a mixed beam of e + and e~ pro- 
duced via pair production from a secondary photon 
beam, was used to simultaneously measure a(e + p— >e + p) 
and cr(e~ p->e~p) elastic scattering cross sections. Cross 
sections in the kinematical range of 0.5 < Q 2 < 2.0 
(GeV/c) 2 and 0.2 < e < 0.9 can be measured. The third 
is the OLYMPUS experiment [H, where the DORIS lep- 
ton storage ring at DESY will be used to extract the ratio 
R e+e - from Q 2 = 0.6 (GeV/c) 2 and e = 0.90, to Q 2 = 
2.2 (GeV/c) 2 and e = 0.35. 

IV. CONCLUSIONS 

In conclusion, we extracted the elastic electromag- 
netic form factors of the proton using two different 
parametrizations for the OPE-TPE interference func- 
tion F(Q 2 ,e): the QA parametrization Eq. (|9]) and 
the BK parametrization Eq. (fTO]) . Both parametriza- 
tions are linear in e, but make different assumptions 
as to where the TPE contributions vanish. In the BK 
parametrization, the TPE correction to <7r was con- 
strained by enforcing the Regge limit which was not 
done in the QA parametrization. In both parametriza- 
tions, we constrained h p Ge p /Gm p using recoil polariza- 
tion data. The values of Gm p and Ge p extracted using 
the QA parametrization are smaller than those obtained 
using the BK parametrization, by as much as 10%, with 
significantly larger uncertainties obtained using the QA 
parametrization. 

The form factor results from the BK fit are generally in 
good agreement with the form factors based on a global 
analysis including calculated TPE corrections [38||, as 
well as some previous phenomenological extractions. For 
Ge p , this is essentially because the polarization trans- 
fer data is assumed to be unaffected by TPE correc- 
tions and taken as a constraint in most of the extrac- 
tions, while for Gm p , the extraction is slightly more 
dependent on the detailed assumptions of the analy- 
ses. Using the BK parametrization, the TPE ampli- 
tude a(Q 2 ) was extracted. The amplitude is on the 
few percent level, and increases in size with increasing 
Q 2 . We parametrize the TPE contribution and fit uncer- 



tainties as a(Q 2 ) = (-0.0191 ± 0.0014) ^/Q 2 ± 0.003 for 
1.5 < Q 2 < 4.0 (GeV/c) 2 . We estimate that the uncer- 
tainties associated with the assumptions made in this ex- 
traction (no TPE correction to polarization data, linear 
of TPE contribution to the cross section, etc.), common 
to most of the extractions presented here, yield an uncer- 
tainty that is comparable to the quoted fit uncertainty, al- 
though these uncertainties are more likely to be strongly 
correlated with Q 2 and of special importance at lower 
Q 2 values. Note that recent high-precision polarization 
measurements of H p Ge p /Gm p at low Q 2 [52H551 ]. com- 
bined with expected results from comparison of positron 
and electron scattering, will significantly improve our 
knowledge of the TPE contributions at lower Q 2 values. 
However, at very low Q 2 values, the TPE contributions 
have significantly more impact on Gm p [56l . [57j . which 
is more difficult to extract precisely in Rosenbluth mea- 
surements. 

Note that the cross sections reported in Ref. [33] were 
determined by detecting recoiling protons, in contrast to 
all other measurements which detected the scattered elec- 
trons. The consistency of the extraction suggests that the 
approximations used to calculate standard radiative cor- 
rections, which yield very different corrections for these 
two cases, are reliable, although the quantitative compar- 
ison is limited by the precision of the electron detection 
measurements. A more extensive set of such measure- 
ments, covering 0.4 < Q 2 <, 5 (GeV/c) 2 is under analysis, 
and will allow for a much more detailed examination [55j . 

Finally, we compare extractions of the TPE contribu- 
tions to the ratio of positron-proton and electron-proton 
scattering cross sections. We use these to make predic- 
tions for the higher Q 2 kinematics of the recently com- 
pleted and ongoing measurements of R e + e -. The lower 
Q 2 values enter the region where the present data are lim- 
ited in their ability to extract TPE contributions based 
on the comparison of cross section and polarization mea- 
surements. The positron measurements will provide new 
information which can provide the first direct experimen- 
tal evidence for TPE contributions at low Q 2 . 
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